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1. Introduction 

A wavelet is a special case of a vector in a separable Hilbert space that generates a 
basis under the action of a collection, or "system" , of unitary operators defined in 
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terms of translation and dilation operations. This approach to wavelet theory goes 
back, in particular, to earlier work of Goodman, Lee and Tang [10] in the context 
of multiresolution analysis. We will begin by describing the operator-interpolation 
approach to wavelet theory using the local commutant of a system that was worked 
out by the speaker and his collaborators a few years ago. This is really an abstract 
application of the theory of operator algebras, mainly von Neumann algebras, 
to wavelet theory. The concrete applications of operator-interpolation to wavelet 
theory include results obtained using specially constructed families of wavelet sets. 
In fact X. Dai and the speaker had originally developed our theory of wavelet sets 
[5] specifically to take advantage of their natural and elegant relationships with 
these wavelet unitary systems. We will also discuss some new results and open 
questions. 

The main idea in operator-theoretic interpolation of wavelets (and frames) 
is that new wavelets can be obtained as linear combinations of known ones using 
coefficients which are operators (in fact, Fourier multipliers) in a certain class. 
Both the ideas and the essential computations extend naturally to more general 
unitary systems and wandering vectors. Many of the methods work for more in- 
volved systems that are important to applied harmonic analysis, such as Gabor 
and generalized Gabor systems, and various types of frame unitary systems. 

1.1. Terminology 

The set of all bounded linear operators on a Hilbert space H will be denoted by 
B{H). A bilateral shift U on H is a unitary operator U for which there exists 
a closed linear subspace E C H with the property that the family of subspaces 
{U"'E : n S Z} are orthogonal and give a direct-sum decomposition of H. The 
subspace E is called a complete wandering subspace for U , and the multiplicity of 
U is defined to be the dimension of E. The strong operator topology on B(H is the 
topology of pointwise convergence, and the weak operator topology is the weakest 
topology such that the vector functionals uJx,y on B{H) defined hy A i-^ {Ax,y), 
A S B{H), x,y € H, are all continuous. An algebra of operators is a linear subspace 
of B{H) which is closed under multiplication. An operator algebra is an algebra 
of operators which is norm-closed. A subset S C B{H) is called selfadjoint if 
whenever A ^ S then also A* ^ S. A C*- algebra is a self-adjoint operator algebra. 
A von Neumann algebra is a C*-algebra which is closed in the weak operator 
topology. For a unital operator algebra, it is well known that being closed in the 
weak operator topology is equivalent to being closed in the closed in the strong 
operator topology. The commutant of a set S of operators in B{H) is the family of 
all operators in B{H) that commute with every operator in S. It is closed under 
addition and multiplication, so is an algebra. And it is clearly closed in both the 
weak operator topology and the strong operator topology. We use the standard 
prime notation for the commutant. So the commutant of a subset S C B{H) is 
denoted: S' := {A G B{H) : AS = SA, S € S}. The commutant of a selfadjoint set 
of operators is clearly a von Neumann algebra. Moreover, by a famous theorem of 
Fuglede every operator which commutes with a normal operator A^ also commutes 
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with its adjoint N* , and hence the comniutant of any set of normal operators is also 
a von Neumann algebra. So, of particular relevance to this work, the comniutant 
of any set of unitary operators is a von Neumann algebra. 

1.1.1. Frames and Operators. A sequence of vectors {fj} in a separable Hilbert 
space iJ is a frame (or frame sequence) if there exist constants Ci , C2 > such 
that 

Ci||/||^<^|</,/,. >p<C2||/||^ 
j 
for all f E H . li Ci — C2 the frame is called tight, and if Ci = C2 = 1, {fj} is 
called a Parseval frame. (The term normalized tight has also been used for this 
(cf [14]). A vector ^ is called a /rame vector for a unitary system U if the set of 
vectors U£, is a frame for H. 

A Riesz basis for a Hilbert space is a bounded unconditional basis. Frames 
sequences are generalizations of Riesz bases. A number of the basic aspects of 
a geometric, or operator-theoretic, approach to discrete frame theory on Hilbert 
space arises from the fact that a frame sequence is simply an "inner" direct sum- 
mand of a Riesz basis. The basic principle is that a Hilbert space frame sequence 
can be dilated to a Riesz basis for larger Hilbert space. We call this the Frame 
Dilation Theorem. In other words, for a given frame sequence there is a larger 
Hilbert space and a Riesz basis for the larger space such that the orthogonal pro- 
jection from the larger space onto the smaller space compresses the Riesz basis to 
the frame sequence. We proved this at the beginning of [14], and used it to prove 
the other results [14], and subsequently to prove some applications to Hilbert C*- 
module theory jointly with M. Frank. We proved it first for Parseval frames, and 
then for general frames. (We remark that this type of dilation result for frames was 
also independently known and used independently by several others in different 
contexts.) 

It is interesting to note that the Parseval frame case of the Frame Dilation 
Theorem can be derived easily from the purely atomic case of a well known theorem 
of Naimark on projection valued measures. We thank Chandler Davis for pointing 
this out to us at the Canadian Operator Algebras Symposium in 1999. We (Han 
and I) basically proved this special case of Naimark's theorem implicitely in the 
first section of [14] without recognizing it was a special case of Naimark's theorem, 
and then we proved the appropriate generalization we needed for general (non- 
tight) frames. Naimark's Dilation Theorem basically states that a suitable positive 
operator valued measure on a Hilbert space dilates to a projection valued measure 
on a larger Hilbert space. That is, there is a projection valued measure [PVM] on a 
larger Hilbert space such that the orthogonal projection from the larger space unto 
the smaller space compresses the [PVM] to the [POVM]. In the discrete (i.e. purely 
atomic measure) case, it can be interpreted as stating that a suitable sequence of 
positive operators dilates to a sequence of projections. The dilation theorem for 
a Parseval frame follows easily from Naimark's Theorem applied to the [POVM] 
obtained by replacing each vector Xi in the frame sequence with the elementary 
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tensor operator Xi®Xi^ obtaining the atoms for a [POVM] defined on all subsets of 
the index set for the frame. The dilation theorem for a general (non-tight) frame 
does not seem to follow directly from Naimark's theorem - but it may follow from 
a generalization of it. We remark that some other generalizations of the frame 
dilation theorem have been recently worked out, notably by W. Czaja. 

1.1.2. Unitary Systems, Wandering Vectors, and Frame Vectors. We define a uni- 
tary system to be simply a countable collection of unitary operators lA acting on 
a Hilbert space H which contains the identity operator. The interesting unitary 
systems all have additional structural properties of various types. (For instance, 
wavelet systems and Gabor systems are both "ordered products" of two abelian 
groups: the dilation and translation groups in the wavelet case, and the modula- 
tion and translation groups in the Gabor case.) We will say that a vector ip Cz H 
is wandering for U if the set 

Uip:^{Uij:U eU} (1) 

is an orthonormal set, and we will call "0 a complete wandering vector for U if 
Uip spans H. This (abstract) point of view can be useful. Write W{U) for the set 
of complete wandering vectors for U. Further, a Riesz vector for hi is a, vector ■0 
such that lA^ is a Riesz basis for H (indexed by the elements oilA), and a frame 
vector is a vector -0 such that lAil) is a frame sequence for H (again using lA as its 
index set), and we adopt similar terminology for Parseval frame vectors and Bessel 
vectors. We use TZW{U), J^{U), VF{U), B{U) to denote, respectively, the sets of 
Riesz vectors, frame vectors, Parseval frame vectors, and Bessel vectors for lA. 

One of the main tools in this work is the local commutant of a system of uni- 
tary operators (see section 3.2). This is a natural generalization of the commutant 
of the system, and like the commutant it is a linear space of operators which is 
closed in the weak and the strong operator topologies, but unlike the commutant it 
is usually not selfadjoint, and is usually not closed under multiplication. It contains 
the commutant of the system, but can be much larger than the commutant. The 
local commutant of a wavelet unitary system captures all the information about 
the wavelet system in an essential way, and this gives the flavor of our approach 
to the subject. 

1.1.3. Normalizers. If J7 is a unitary operator and A is an operator algebra, then 
U is said to normalize A if U* -A-U = A .In the most interesting cases of operator- 
theoretic interpolation: that is, for those cases that yield the strongest structural 
results, the relevant unitaries in the local commutant of the system normalize the 
commutant of the system. 

1.2. Acknowledgement 

This article was written in response to an invitation by the organizers of the 4th 
International Conference on Wavelet Analysis and its Applications, Macau, China, 
December 2005 [WAA2005] to be a keynote speaker. We thank the organizers for 
their kind invitation to present these notes. 
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2. Wavelets 



For simplicity of presentation, much of the work in this article will deal with one- 
dimensional wavelets, and in particular, the dyadic case. The other cases: non- 
dyadic wavelets and wavelets in higher dimensions, are at least notationally more 
complicated. 

2.1. One Dimension 

A dyadic orthonormal wavelet in one dimension is a unit vector ijj G Lp'iM.,^), 
with /i Lebesgue measure, with the property that the set 

{2^V'(2"i-0 :",'e Z} (2) 

of all integral translates of '0 followed by dilations by arbitrary integral powers 
of 2, is an orthonormal basis for L^(E,/i). The term dyadic refers to the dilation 
factor "2". The term mother wavelet is also used in the literature for -0. Then the 
functions 

xPnj :=2t0(2"t-O 
are called elements of the wavelet basis generated by the " mother" . The functions 
"071,; will not themselves be mother wavelets unless n — Q. 

Let T and D be the translation (by 1) and dilation (by 2) unitary operators 
in B(i2(K) given by {Tf)(t) = f{t - 1) and {Df){t) = V2f{2t). Then 

25V(2"t-0 = {D'^T^%l:){t) 

for all n, / e Z. Operator-theoretically, the operators T, D are bilateral shifts of 
infinite multiplicity. It is obvious that i^([0, 1]), considered as a subspace of L^(IR), 
is a complete wandering subspace for T, and that L^([— 2, — 1] U [1, 2]) is a complete 
wandering subspace for D. 

Let hiD,T be the unitary system defined by 

Ud,t = {£'"r' : n, / e Z} (3) 

where D and T are the operators defined above. Then ^/^ is a dyadic orthonormal 
wavelet if and only if -0 is a complete wandering vector for the unitary system 
l^D,T- This was our original motivation for developing the abstract unitary system 
theory. Write 

W{D,T):^W{Ud,t) (4) 

to denote the set of all dyadic orthonormal wavelets in one dimension. 

An abstract interpretation is that, since £> is a bilateral shift it has (many) 
complete wandering subspaces, and a wavelet for the system is a vector -0 whose 
translation space (that is, the closed linear span of {T^ : fc G Z} is a complete 
wandering subspace for D. Hence V' must generate an orthonormal basis for the 
entire Hilbert space under the action of the unitary system. 

In one dimension, there are non-dyadic orthonormal wavelets: i.e. wavelets 
for all possible dilation factors besides 2 (the dyadic case). We said "possible", 
because the scales {0, 1, —1} are excluded as scales because the dilation operators 
they would introduce are not bilateral shifts. All other real numbers for scales 



6 Larson 

yield wavelet theories. In [5, Example 4.5 (x)] a family of examples is given of 
three-interval wavelet sets (and hence wavelets) for all scales d > 2, and it was 
noted there that such a family also exists for dilation factors 1 < d < 2. There 
is some recent (yet unpublished) work that has been done, by REU students and 
mentors, building on this, classifying finite-interval wavelet sets for all possible real 
(positive and negative scale factors). I will mention this work, in passing, in my 
talk. 

2.2. N-Dimensions 

2.2.1. Expansive Dilations. Let 1 < m < oo, and let ^ be an n x n real matrix 
which is expansive (equivalently, all (complex) eigenvalues have modulus > 1). 
By a dilation - A regular-translation orthonormal wavelet we mean a function 
V- e L2(M") such that 

{\detiApi;{A^t - (WM. -, Inf) -.n^leZ} (5) 

where t = {ti, ...,t„)*, is an orthonormal basis for i^(M";r7i). (Here m is product 
Lebesgue measure, and the superscript "t" means transpose.) 

li A € M„(]R) is invertible (so in particular if A is expansive), then it is very 
easy to verify that the operator defined by 

iDAf)(t) = \detA\if{At) (6) 

for / G L^(M"), t G K", is unitary. For 1 < i < n, let Ti be the unitary operator 
determined by translation by 1 in the i*^ coordinate direction. The set (5) above 
is then 

{D'XTi^---T^-i^:k,heZ} (7) 

If the dilation matrix A is expansive, but the translations are along some 
oblique lattice, then there is an invertible real nxn matrix T such that conjugation 
with Dt takes the entire wavelet system to a regular-translation expansive-dilation 
matrix. This is easily worked out, and was shown in detail in [18] in the context of 
working out a complete theory of unitary equivalence of wavelet systems. Hence 
the wavelet theories are equivalent. 

2.2.2. Non-Expansive Dilations. Much work has been accomplished concerning 
the existence of wavelets for dilation matrices A which are not expansive. Some of 
the original work was accomplished in the Ph.D. theses of Q. Gu and D. Speegle, 
when they were together finishing up at Texas A&M. Some significant additional 
work was accomplished by Speegle and also by others. In [18], with lonascu and 
Pearcy we proved that if an nxn real invertible matrix A is not similar (in the nxn 
complex matrices) to a unitary matrix, then the corresponding dilation operator 
Da is in fact a bilateral shift of infinite multiplicity. If a dilation matrix were to 
admit any type of wavelet (or frame- wavelet) theory, then it is well-known that a 
necessary condition would be that the corresponding dilation operator would have 
to be a bilateral shift of infinite multiplicity. I am happy to report that in very 
recent work [23], with E. Schulz, D. Speegle, and K. Taylor, we have succeeded 
in showing that this minimal condition is in fact sufficient: such a matrix, with 
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regular translation lattice, admits a (perhaps infinite) tuple of functions, which 
collectively generates a frame-wavelet under the action of this unitary system. 

3. More General Unitary Systems 

3.1. Some Restrictions 

We note that most unitary systems U do not have complete wandering vectors. 
For W{U) to be nonempty, the set U must be very special. It must be countable if 
it acts separably (i.e. on a separable Hilbert space), and it must be discrete in the 
strong operator topology because ii U,V Cz U and if x is a wandering vector for U 
then 

\\U-V\\ > \\Ux-Vx\\ ^V2 

Certain other properties are forced on U by the presence of a wandering vector. 
(Or indeed, by the nontriviality of any of the sets W(U), nW{U), T{U), VT{U), 
B{U).) One purpose of [5] was to investigate such properties. Indeed, it was a 
matter of some surprise to us to discover that such a theory is viable even in some 
considerable generality. For perspective, it is useful to note that while Ud.t has 
complete wandering vectors, the reversed system 

Kt^d == {T'D" : n, / e Z} 

fails to have a complete wandering vector. (A proof of this was given in the intro- 
duction to [5].) 

3.2. The Local Commutant 

3.2.1. A Special Case: The System Ud,t- Computational aspects of operator the- 
ory can be introduced into the wavelet framework in an elementary way. Here is 
the way we originally did it: Fix a wavelet ip and consider the set of all operators 
A e _B(L^(]R)) which commute with the action of dilation and translation on %Ij. 
That is, require 

(AV')(2"f-0 = A(^(2"t-0) (8) 

or equivalently 

D'^T^A^ ^ AD'^T^tP (9) 

for all n, / G Z. Call this the local commutant of the wavelet system Ud.t o-t 
the vector -0. (In our first preliminary writings and talks we called it the point 
commutant of the system.) Formally, the local commutant of the dyadic wavelet 
system on L^(R) is: 

C-4,{Ud,t) :={AeB(L2(R)) : (AL>"T' - L>"r'A)?A = 0, Vn, ? e Z} (10) 

This is a linear subspace of B{H) which is closed in the strong operator topol- 
ogy, and in the weak operator topology, and it clearly contains the commutant of 
{AT}. 

A motivating example is that if 77 is any other wavelet, let V :— Vjj be the 
unitary (we call it the interpolation unitary) that takes the basis '>pn,i to the basis 
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r]n,i. That is, V^pnj = Vnd for all n, ? e Z. Then 77 = V^j, so VD'^T^iP = D'^T^Vi: 
hence F G C^{Ud,t)- 

In the case of a pair of complete wandering vectors if), rj for a general unitary 
system lA, we will use the same notation VV for the unitary that takes the vector 
Ui) to Uri for all U eU. 

This simple-minded idea is reversible, so for every unitary V in C^{Ud.t) 
the vector V-ip is a wavelet. This correspondence between unitaries in C^{D,T) 
and dyadic orthonormal wavelets is one-to-one and onto (see Proposition 3.1). 
This turns out to be useful, because it leads to some new formulas relating to 
decomposition and factorization results for wavelets, making use of the linear and 
multiplicative properties oi C^{D,T). 

It turns out (a proof is required) that the entire local commutant of the sys- 
tem IAd,t at a wavelet V' is not closed under multiplication, but it also turns out 
(also via a proof) that for most (and perhaps all) wavelets ■!/; the local commutant 
at "0 contains many noncommutative operator algebras (in fact von Neumann al- 
gebras) as subsets, and their unitary groups parameterize norm-arcwise-connected 
families of wavelets. Moreover, C^{D,T) is closed under left m,ultiplication by the 
commutant {D, T}', which turns out to be an abelian nonatomic von Neumann 
algebra. The fact that C^{D,T) is a left module under {D,T}' leads to a method 
of obtaining new wavelets from old, and of obtaining connectedness results for 
wavelets, which we called operator-theoretic interpolation of wavelets in [5], (or 
simply operator-interpolation). 

3.2.2. General Systems. More generally, let S C B{H) be a set of operators, where 
H \s a. separable Hilbert space, and let x G iJ be a nonzero vector, and formally 
define the local commutant of 5 at x by 

C^{S) := {A e B{H) : {AS - SA)x = 0, 5 G 5} 

As in the wavelet case, this is a weakly and strongly closed linear subspace of 
B{H) which contains the commutant S' of S. If x is cyclic for S in the sense that 
span(5a;) is dense in H, then x separates Cx{S) in the sense that for S G Cx{S), 
we have Sx = iff x = 0. Indeed, if A G Cx{S) and if Ax = 0, then for any S G S 
we have ASx = SAx = 0, so ASx = 0, and hence A = 0. 

If A G Cx{S) and B G 5', let C = BA. Then for aU 5 G 5, 

(CS - SC)x = B{AS)x - {SB)Ax ^ B{SA)x - {BS)Ax = 

because ASx — SAx since A G Cx{S), and SB — BS since B G S' . Hence Cx{S) is 
closed under left multiplication by operators in S'. That is, Cx{S is a left module 
over S' . 

It is interesting that, if in addition 5 is a multiplicative semigroup, then in 
fact Cx{S) is identical with the commutant S' so in this case the commutant is not 
a new structure. To see this, suppose A G Cx{S). Then for each S,T G S we have 
ST G S, and so 

AS{Tx) == {ST)Ax ^ S{ATx) ^ {S)Tx 
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So since T Cz S was arbitrary and span(5a;) = H, it follows that AS = SA. 

Proposition 3.1. If U is any unitary system for which W{U) ^ 0, then for any 
tp e W{U) 

WilA) = {Uip : U is a unitary operator in C,p(U)} 
and the correspondence U -^ Uip is one-to-one. 

Proposition 3.2. LetlA be a unitary system on a Hilbert space H . If-tp is a complete 
wandering vector for hi , then: 

(i) TZWihi) — {Alp : A is an operator in C^ilA) that is invertible in B{H)}; 

(ii) T{U) = {Alp : A is an operator in C^(U) that is surjective}; 

(ii) VT{IA) = {Alp : A is an operator in C^(U) that is a co-isometry}; 

(ii) B{U) = {Alp : A is an operator in C^{U)} 

3.3. Operator-Theoretic Interpolation 

Now suppose U is a unitary system, such as Ud,t, and suppose {ipi, tp2, ■ ■ ■ , V'm} C 
W{U). (In the case oi Ud,t, this means that (ipi,tp2, ■ ■ ■ jipn) is an n-tuple of 
wavelets. 

Let (Ai, A2, . . . , An) be an n-tuple of operators in the commutant W of U, 
and let rj be the vector 

Then 

= [A, + A2V^^ + ■■■ + A„T/;;)Vi . (11) 

We say that rj is obtained by operator interpolation from {tpi, 4'2, ■ ■ ■ , V'm}- 
Since C^^ iU) is a left U' - module, it follows that the operator 

A:^A,+A2V^-^+...AnV^- (12) 

is an clement of C^^iU). Moreover, if B is another element of C^^iU) such that 
f] = Btpi, then A — B e C^^ {U) and [A — B)^pi — 0. So since ipi separates C^^ (U) 
it follows that A = B. Thus A is the unique element oiC^-^{U) that takes ipi to rj. 
Let 5^i,....^„ be the family of all finite sums of the form 

n 

This is the left module of W generated by {/, Vf^, . . . , V^^}. It is the W -linear 

spanof{7,y;;, •■•,<;}• 
Let 

A^l/'l,•••.'/'-^ •= ('5^1,...,^„)V'l (13) 
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So 

We call this the interpolation space for U generated by (-01, . . . , ^„). From the 
above discussion, it follows that for every vector r) G ■A4ipi,^2,---,tl'n there exists a 
unique operator A G C^^ {U) such that rj = Aipi, and moreover this A is an element 

of '5^i,...,i/.„- 

3.3.1. Normalizing the Commutant. In certain essential cases (and we are not sure 
how general this type of case is) one can prove that an interpolation unitary V^ 
normalizes the commutant W of the system in the sense that V^U'V2 = W . (Here, 
it is easily seen that {V^)* = V^ .) Write V := V^- If V normalizes W , then the 
algebra, before norm closure, generated by hi' and V is the set of all finite sums 
(trig polynomials) of the form ^ A„y", with coefficients An G W , n G Z. The 
closure in the strong operator topology is a von Neumann algebra. Now suppose 
further that every power of V is contained in C^{IA). This occurs only in special 
cases, yet it occurs frequently enough to yield some general methods. Then since 
C^{U) is a SOT-closed linear subspace which is closed under left multiplication 
by W , this von Neumann algebra is contained in C^{IA)^ so its unitary group 
parameterizes a norm-path-connected subset of W{IA) that contains -0 and rj via 
the correspondence U — > Ui}}. 

In the special case of wavelets, this is the basis for the work that Dai and I 
did in [5, Chapter 5] on operator-theoretic interpolation of wavelets. In fact, we 
specialized there and reserved the term operator-theoretic interpolation to refer 
explicitly to the case when the interpolation unitaries normalize the commutant. 
In some subsequent work, we loosened this restriction yielding our more general 
definition given in this article, because there are cases of interest in which we 
weren't able to prove normalization. However, it turns out that if -0 and rj are 
s-elementary wavelets (see section 4.4), then indeed V^ normalizes {Z?,r}'. (See 
Proposition 5.3.) Moreover, V^ has a very special form: after conjugating with 
the Fourier transform, it is a composition operator with a symbol a that is a 
natural and very computable measure-preserving transformation of R. In fact, it 
is precisely this special form for V2 that allows us to make the computation that 
it normalizes {D, T}' . On the other hand, we know of no pair (-0, rf) of wavelets for 
which V2 fails to normalize {£),T}'. The difficulty is simply that in general it is 
very hard to do the computations. This is stated as Problem 2 in the final section 
on Open Problems. 

In the wavelet case Ud.t , if "0 G W{D, T) then it turns out that C^{Ud.t) 
is in fact much larger than {Ud,t)' — {D, T}' , underscoring the fact that Ud,t is 
NOT a group. In particular, {D,r}' is abelian while C^{'D,T) is nonabelian for 
every wavelet ip. (The proof of these facts are contained in [5].) 
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3.3.2. Interpolation Pairs of Wandering Vectors. In some cases where a pair V', V 
of vectors in yV{U) are given it turns out that the unitary V in C^{U) with Vip = 
rj happens to be a symmetry (i.e. V'^ = I). Such pairs are caUed interpolation 
pairs of wandering vectors, and in the case where U is a wavelet system, they 
are called interpolation pairs of wavelets. Interpolation pairs are more prevalent 
in the theory, and in particular the wavelet theory, than one might expect. In 
this case (and in more complex generalizations of this) certain linear combinations 
of complete wandering vectors are themselves complete wandering vectors ~ not 
simply complete Riesz vectors. 

Proposition 3.3. Let U be a unitary system, let ip,ri £ W{U), and let V be the 
unique operator in C^{U) with Vip = rj. Suppose 

1/2 = /. 

Then 

cos a ■ -ip + i sin a ■ rj £ W{U) 
for alio <a< lix . 

The above result can be thought of as the prototype of our operator-theoretic 
interpolation results. It is the second most elementary case. (The most elementary 
case is described in the context of the exposition of Problem 4 in the final section.) 
More generally, the scalar a in Proposition 3.3 can be replaced with an appropriate 
self-adjoint operator in the commutant oiU. In the wavelet case, after conjugating 
with the Fourier transform, which is a unitary operator, this means that a can be 
replaced with a wide class of nonnegative dilation-periodic (see definition below) 
bounded measurable functions on R. 



4. Wavelet Sets 

Wavelet sets belong to the theory of wavelets via the Fourier Transform. We will 
do most of this section in a tutorial-style, to make the concepts more accessible to 
students and colleagues who are not already familiar with them. 

4.1. Fourier Transform 

We will use the following form of the Fourier Plancherel transform T on Ti = 
i2(R), which is a form that is normalized so it is a unitary transformation, a 
property that is desirable for our treatment. 
If /,ge Li(R)nL2(M) then 

{Tf){s) := -^ / e'^''f{t)dt := />), (14) 



/27r 
and 



{T-'g){t) = -= / e^^'g{s)ds. (15) 

y Ait 
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We have 

V ^TT JK 

So TTaT-^g = e-*''"^. For A e B{n) let i denote J^AJ^-^. Thus 

f„=Me-, = a, (16) 

where for h E L°^ we use M/j to denote the multiphcation operator f ^ hf. Since 
{Me-.== : a e M} generates the m.a.s.a. X'(M) := {Mh: h e L°°(R)} as a von 
Neumann algebra, we have 

J'AtJ''^ =V{R). 
Similarly, 

{TD'\f){s) = -L / e-^^*(\/2)"/(2"t)di 



'2tt 

= (a/2)-" . -^ / e-^'-''^'f{t)dt 

V 27r jR 

= (V2)-2(^/)(22-"^) = (i?-"^/)(s). 

So 5" = £>-" = D™. Therefore, 

D = D-^=D*. (17) 

If / is an _L^(]R) function, as usual we write f{s) = {T{f)){s). If p{s) is a 
real-valued function such that f{s) — e*'''^''^|/(s)|, we call p{s) the phase of/. The 
phase is well defined a.e. modulo 27r-translation. 

4.2. The Commutant of {D,T} 

We have T{D,TyT-'^ = {S,f}'. It turns out that {D,f}' has an elementary 
characterization in terms of Fourier multipliers: 

Theorem 4.1. 

{D,fy ^ {Mh-. h e L°°(M) and h{s) ^ h{2s) a.e.}. 

Proof. Since D — D* and D is unitary, it is clear that Mh E {D, T}' if and only 
if Mh commutes with D. So let g E i^(M) be arbitrary. Then (a.e.) we have 

iMhDg){s) = h{s)(V2 g{2s)), and 

iDMhg){s) = D{h{s)g{s)) - V7I(2s)g(2s). 

Since these must be equal a.e. for arbitrary g, we must have h{s) = h{2s) a.e. D 

Now let E = [-2, -1)U [1, 2), and forn e Z let E„ = {2"a;: x E E}. Observe 
that the sets E„ are disjoint and have union R\{0}. So if g is any uniformly 
bounded function on E, then g extends uniquely (a.e.) to a function g E L°°{M.) 
satisfying 

g{s) = g{2s), s E K, 
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by setting 

g(2"s) := 5(s), seE,neZ, 

and g{0) — 0. We have \\g\\oo = Halloo- Conversely, if h is any function satisfying 
h{s) — h(2s) a.e., then h is uniquely (a.e.) determined by its restriction to E. This 
1-1 mapping g -^ Mg from L°°{E) onto {D^T}' is a ^-isomorphism. 

We will refer to a function h satisfying h{s) = h{2s) a.e. as a 2- dilation 
periodic function. This gives a simple algorithm for computing a large class of 
wavelets from a given one, by simply modifying the phase (see also section 4.7): 

Given ip,, let ip — J^{ip), choose a real-valued function h G L°°{E) 
arbitrarily, let g — exp(ih), extend to a 2-dilation periodic (18) 

function g as above, and compute tpg — T^ (ffV')- 

In the description above, the set E could clearly be replaced with [— 27r, — 7r)U 
[tt, 27r), or with any other "dyadic" set [—2a, a) U [a, 2a) for some a > 0. 

4.3. The Shannon Wavelet 

We now give an account of s-elementary and MS* i^- wavelets. The two most ele- 
mentary dyadic orthonormal wavelets are the well-known Haar wavelet and Shan- 
non's wavelet (also called the Littlewood-Paley wavelet). The Shannon set is the 
prototype of the class of wavelet sets. 

Shannon's wavelet is the L^(IR)-function with Fourier transform ^5 — -y=XEo 
where 

£^0 = [-27r,-7r)U [7r,27r). (19) 

The argument that -05 is a wavelet is in a way even more transparent than for 
the Haar wavelet. And it has the advantage of generalizing nicely. For a simple 
argument, start from the fact that the exponents 

restricted to [0, 27r] and normalized by ^= is an orthonormal basis for £^[0, 27r]. 



e 



Write £^0 = £-U£+ where £_ = [-27r, -tt), £+ = [7r,27r). Since {£_ +27r, £+} is 
a partition of [0, 27r) and since the exponentials e*^^ are invariant under translation 
by 27r, it follows that 

„its 

neZ} (20) 



/27r 
is an orthonormal basis for L'^(Eo). Since T — M^-is, this set can be written 

{f%: £eZ}. (21) 

Next, note that any "dyadic interval" of the form J = [b, 26), for some 6 > has 
the property that {2" J: n € Z}, is a partition of (0, 00). Similarly, any set of the 
form 

JC=[-2a,-a)U[b,2b) (22) 
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for a, 6 > 0, has the property that 

{2"/C: neZ} 

is a partition of ]R\{0}. It follows that the space L'^(JC), considered as a subspace 
of L^(R), is a complete wandering subspace for the dilation unitary (Df){s) ~ 
y/2 /(2s). For each n G Z, 

^"(^^(/C)) =l2(2-"/C). (23) 

So ®„ -D"(L^(/C)) is a direct sum decomposition of L^(K). In particular Eq has 
this property. So 

eez} (24) 



D" 



£ez: ■ ^ 



2-iBo 



So J L vStt 

is an orthonormal basis for L^(2^"_Eo) for each n. It follows that 

is an orthonormal basis for L^(M). Hence {D^T^ipg ■ n,^ G Z} is an orthonormal 
basis for L^(M), as required. 

For our work, in order to proceed with developing an operator-algebraic the- 
ory that had a chance of directly impacting concrete function-theoretic wavelet 
theory we needed a large supply of examples of wavelets which were elementary 
enough to work with. First, we found another "Shannon- type" wavelet in the 
literature. This was the Journe wavelet, which we found described on p. 136 in 
Daubechies book [8]. Its Fourier transform is 'ipj = -^=Xe,j, where 



Ej = 



327r ^ 

,-47r 

7 ' 



_ 1 

V2tt 


XEj 


"47r \ 


U 





47r\ 4tt \ , 327r 

U — ,7r U 47r, 

7 7 ' '7 



Then, thinking the old adage "where there's smoke there's fire!" , we painstakingly 
worked out many more examples. So far, these are the basic building blocks in the 
concrete part of our theory. By this we mean the part of our theory that has had 
some type of direct impact on function-theoretic wavelet theory. 

4.4. Definition of Wavelet Set 

We define a wavelet set to be a measurable subset i? of M for which — ^vc; is the 
Fourier transform of a wavelet. The wavelet V'e '■— -^=Xe is called s- elementary 
in [5]. 

It turns out that this class of wavelets was also discovered and systematically 
explored completely independently, and in about the same time period, by Guido 
Weiss (Washington University), his colleague and former student E. Hernandez (U. 
Madrid), and his students X. Fang and X. Wang. Two of the papers of this group 
are [9] and [17], in which they are called MSF (minimally supported frequency) 
wavelets. In signal processing, the parameter s, which is the independent vari- 
able for ip, is the frequency variable, and the variable t, which is the independent 
variable for ■!/;, is the time variable. It is not hard to show that no function with 
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support a subset of a wavelet set E of strictly smaller measure can be the Fourier 
transform of a wavelet. (Here, the support of a measurable function is defined to 
be the set of points at which it does not vanish.) In other words, an MSF wavelet 
has minimal possible support in the frequency domain. However, the problem of 
whether the support set of any wavelet necessarily contains a wavelet set remains 
open. It was raised by this author (Larson) in a talk about ten years ago, and has 
been open for several years. We include it as Problem 3 in the final section of this 
article. A natural subproblem, which was posed in the same talk, asks whether a 
wavelet with minimal possible support in the frequency domain is in fact an MSF 
wavelet; or equivalently, is its support a wavelet set? 

4.4.1. The Spectral Set Condition. ^From the argument above describing why 
Shannon's wavelet is, indeed, a wavelet, it is clear that sufficient conditions for E 
to be a wavelet set are 

(i) the normalized exponential —h=e^^'^^ ^ € Z, when restricted to E 

should constitute an orthonormal basis for L'^{E) (in other words E is 
a spectral set for the integer lattice Z), 

and 

(ii) The family {2"E: n G Z} of dilates of E by integral powers of 2 
should constitute a measurable partition (i.e. a partition modulo null 
sets) of R. 

These conditions are also necessary. In fact if a set E satisfies (i), then for it to 
be a wavelet set it is obvious that (ii) must be satisfied. To show that (i) must be 
satisfied by a wavelet set E, consider the vectors 

^""^^ " ^^X2-"B, neZ. 



Since tpE is a wavelet these must be orthogonal, and so the sets {2"E: n E Z} 
must be disjoint modulo null sets. It follows that {—^c'^^'^Ie'- i? G Z} is not only 
an orthonormal set of vectors in L'^{E), it must also span Lp'{E). 

It is known from the theory of spectral sets (as an elementary special case) 
that a measurable set E satisfies (i) if and only if it is a generator of a measurable 
partition of M under translation by 27r (i.e. iff {E + 27rn : n € Z} is a measurable 
partition of R) . This result generalizes to spectral sets for the integral lattice in 
K". For this elementary special case a direct proof is not hard. 

4.5. Translation and Dilation Congruence 

We say that measurable sets E, F are translation congruent modulo lix if there is 
a measurable bijection 0: E -^ F such that (/)(s) — s is an integral multiple of 27r 
for each s d E; ot equivalently, if there is a measurable partition {En '■ n € Z} of 
E such that 

{En + 2mT: n G Z} (25) 
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is a measurable partition of F. Analogously, define measurable sets G and H to 
be dilation congruent modulo 2 if there is a measurable bijection t : G ^ H such 
that for each s G G there is an integer n, depending on s, such that t{s) — 2"s; 
or equivalently, if there is a measurable partition {Gn}?foo of ^ such that 

{2"G}-^ (26) 

is a measurable partition of H. (Translation and dilation congruency modulo other 
positive numbers of course make sense as well.) 
The following lemma is useful. 

Lemma 4.2. Let f G L^(R), and let E = supp(/). Then f has the property that 

is an orthonornial basis for L^ (E) if and only if 
(i) E is congruent to [0, 27r) modulo 27r, and 
(ii) \.fi^)\^'-h «-e- onE. 



27r 



If i<^ is a measurable set which is 27r-translation congruent to [0, 27r), then 
since 



pi£s 



[0,27r) 

is an orthonornial basis for L^[0, 2tt] and the exponentials e are 27r- invariant, as 
in the case of Shannon's wavelet it follows that 



gj£s 



: £e 

E 

is an orthonormal basis for L?{E). Also, if i? is 27r-translation congruent to [0, 27r), 
then since 

{[0,27r) + 27m: n G Z} 

is a measurable partition of R, so is 

{E + 2TTn: nG Z}. 

These arguments can be reversed. 

We say that a measurable subset G C M. is a 2- dilation generator of a partition 
of M if the sets 

2"G:={2"s: s G G}, nGZ (27) 

are disjoint and M\U„2"G is a null set. Also, we say that i? C K is a 2-11 -translation 
generator of a partition of M if the sets 

E + 2n'K := {s + 2m: : seE), nGZ, (28) 

are disjoint and K\ U„ {E + 2n7r) is a null set. 
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Lemma 4.3. A measurable set E Q R is a 2t: -translation generator of a partition 
ofM. if and only if, modulo a null set, E is translation congruent to [0, 27r) modulo 
2iT. Also, a measurable set G C M is a 2-dilation generator of a partition of M 
if and only if, modulo a null set, G is a dilation congruent modulo 2 to the set 

[-27r, -7r)U [tt, 27r). 

4.6. A Criterion 

The following is a useful criterion for wavelet sets. It was published independently 
by Dai-Larson in [5] and by Fang- Wang (who were students of Guido Weiss) in 
[9] at about the same time (in December, 1994). In fact, it is amusing that the two 
papers had been submitted within two days of each other; only much later did we 
even learn of each other's work on wavelets and of this incredible timing. 

Proposition 4.4. Let E <Z M. be a measurable set. Then E is a wavelet set if and 
only if E is both a 2-dilation generator of a partition (modulo null sets) of R and 
a 2t: -translation generator of a partition (modulo null sets) o/R. Equivalently, E 
is a wavelet set if and only if E is both translation congruent to [0, 27r) modulo 27r 
and dilation congruent to [— 27r, —vr) U [tt, 27r) modulo 2. 

Note that a set is 27r-translation congruent to [0, 27r) iff it is 27r-translation 
congruent to [— 27r, tt) U [tt, 27r). So the last sentence of Proposition 14.41 can be 
stated: A measurable set £■ is a wavelet set if and only if it is both 27r-translation 
and 2-dilation congruent to the Littlewood-Paley set [— 27r, — tt) U [tt, 2tt). 

4.7. Phases 

If _E is a wavelet set, and if f{s) is any function with support E which has constant 
modulus ^= on E, then T^^{f) is a wavelet. Indeed, by Lemma H^ IT^ f : G Z} 
is an orthonormal basis for L'^{E), and since the sets 2"£' partition R, so L'^{E) 
is a complete wandering subspace for D, it follows that {D"T^f: n,£ E Z} must 
be an orthonormal basis for L^(R), as required. In [9, 17] the term MSF-wavelet 
includes this type of wavelet. So MSF-wavelets can have arbitrary phase and s- 
elementary wavelets have phase 0. If V" is a wavelet we say [5] that a real- valued 
function p(s) is attainable as a phase of ijj if the function e'''^^''|^/'(s)| is also the 
Fourier transform of a wavelet. So every phase is attainable in this sense for an 
MSF or s-elementary wavelet. Attainable phases of wavelets have been studied in 
[5] and [26], in particular. 

4.8. Some Examples of One-Dimensional Wavelet Sets 

It is usually easy to determine, using the dilation-translation criteria, in Propo- 
sition ^31 whether a given finite union of intervals is a wavelet set. In fact, to 
verify that a given "candidate" set E is a. wavelet set, it is clear from the above 
discussion and criteria that it suffices to do two things. 

(1) Show, by appropriate partitioning, that E is 2-dilation-congruent to 

a set of the form [—2a, —a) U [b, 26) for some a,b > 0. 
and 
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(2) Show, by appropriate partitioning, that E is 27r-translation-congruent 
to a set of the form [c, c + 2tt) for sonic real number c. 

On the other hand, wavelet sets suitable for testing hypotheses can be quite 
difficult to construct. There are very few "recipes" for wavelet sets, as it were. 
Many families of such sets have been constructed for reasons including perspective, 
experimentation, testing hypotheses, etc., including perhaps the pure enjoyment of 
doing the computations — which are somewhat "puzzle-like" in nature. In working 
with the theory it is nice (and in fact we find it necessary) to have a large supply 
of wavelets on hand that permit relatively simple analysis. 

For this reason we take the opportunity here to present for the reader a 
collection of such sets, mainly taken from [5], leaving most of the work in verifying 
that they are indeed wavelet sets to the reader. 

We refer the reader to [6] for a proof of the existence of wavelet sets in K^") , 
and a proof that there are sufficiently many to generate the Borel structure of 
]]j(n) These results are true for arbitrary expansive dilation factors. Some concrete 
examples in the plane were subsequently obtained by Soardi and Weiland, and 
others were obtained by Gu and by Speegle in their thesis work at A&M. Two had 
also been obtained by Dai for inclusion in the revised concluding remarks section 
of our Memoir [5] . 

In these examples we will usually write intervals as half-open intervals [ • , ) 
because it is easier to verify the translation and dilation congruency relations (1) 
and (2) above when wavelet sets are written thus, even though in actuality the 
relations need only hold modulo null sets. 

(i) As mentioned above, an example due to Journe of a wavelet which admits 
no multiresolution analysis is the s-elementary wavelet with wavelet set 



327r 



-4tt U 



47r 
7 



-TT,— U 



47r 
7 



^,7r U 



. 327r 

47r, 

' 7 



327r 


167r^ 


U 


■IGtt 327r 


7 ' 


7 J 


7 ' 7 



To see that this satisfies the criteria, label these intervals, in order, as Ji, J2, J3, Ji 
and write J — UJ,. Then 



Ji U 4 J2 U 4 J3 U J4 = 



This has the form [—2a, a) U [h, 26) so is a 2-dilation generator of a partition of 
M\{0}. Then also observe that 

{ Ji + 67r, J2 + 27r, J3, J4 - 47r} 

is a partition of [0, 27r). 

(ii) The Shannon (or Littlcwood-Paley) set can be generalized. For any — tt < 
a < TT, the set 

Ea = [-27r + 2a, -tt + a) U [tt + a, 27r + 2a) 

is a wavelet set. Indeed, it is clearly a 2-dilation generator of a partition of K\{0}, 
and to see that it satisfies the translation congruency criterion for — tt < a < Q 
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(the case < a < tt is analogous) just observe that 

{[-2tt + 2a, 27r) + 4tt, [-27r, -n + a) + 27r, [tt + a, 27r + 2a)} 

is a partition of [0, 27r). It is clear that ipE^ is then a continuous (in i^(M)-norni) 
path of s-elementary wavelets. Note that 



lim -tpE^ 



/2n 



X[2ttAtt)- 



This is not the Fourier transform of a wavelet because the set [27r, 47r) is not a 
2-dilation generator of a partition of ]R\{0}. So 

lim tpEa 

is not an orthogonal wavelet. (It is what is known as a Hardy wavelet because it 
generates an orthonormal basis for H'^{M.) under dilation and translation.) This 
example demonstrates that W(D,T) is not closed in i^(K). 

(iii) Journe's example above can be extended to a path. For — y < /? < y the 
set 



J6 = 



327T 

^,-47r + 4/3)U 



An , 
-n + (3,-—]U 



47r 

y,7r + /?) U 



47r 

47r + 4/3, 47r H 

7 



is a wavelet set. The same argument in (i) establishes dilation congruency. For 
translation, the argument in (i) shows congruency to [4/3, 27r + 4/3) which is in turn 
congruent to [0, 27r) as required. Observe that here, as opposed to in (ii) above, 
the limit of 'ipj„ as /3 approaches the boundary point j is a wavelet. Its wavelet 
set is a union of 3 disjoint intervals. 

(iv) Let A C [tt, ^) be an arbitrary measurable subset. Then there is a 



wavelet set W, such that W n [tt, ^] 



A. For the construction, let 



B = [27r, 37r)\2A, 



C 



-n,-^)\{A-27r) 



and D = 2A- An. 



Let 



W ^ 



— ,2tt 1 UAUBUCUD. 



We have W n [tt, ^) = A. Observe that the sets [^, 27r), A, B, C, D, arc disjoint. 
Also observe that the sets 

^,2ttYa,^B,2C,D, 

are disjoint and have union [— 27r, — tt) U [tt, 27r). In addition, observe that the sets 

"Stt \ 
— , 27r j , ^, B - 27r, C + 27r, D + 27r, 

are disjoint and have union [0, 27r). Hence VF is a wavelet set. 
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(v) Wavelet sets for arbitrary (not necessarily integral) dilation factors other 
then 2 exist. For instance, if d > 2 is arbitrary, let 

2dTT 2tt 



A = 
B ^ 



d+V d+1 

27r 2tt 



d^ -r d + 1 

2dTT 2d'^iT 



1 ' rf2 _ 1 



and let G = A\JB\JC. Then G is d-wavelet set. To see this, note that {A+2t:, B, G} 
is a partition of an interval of length 27r. So G is 27r-translation-congruent to [0, 27r). 
Also, {A, B, d~^G} is a partition of the set [~da, —a) U [/?, d/3) for a = -^p^, and 
j3 = ^jTrj, so from this form it follows that {d^'^G: n <E Z} is a partition of 
R\{0}. Hence if V' := -^"^(tti^Xg), it follows that {d'^i;{d''t - t): n,€ S Z} is 
orthonormal basis for Lp'i^), as required. 



5. Operator- Theoretic Interpolation for Wavelet Sets 

Operator-theoretic interpolation takes a particularly natural form for the special 
case of s-elementary (or MSF) wavelets that facilitates hands-on computational 
techniques in investigating its properties. Let E^F he a pair of wavelet sets. Then 
for (a.e.) x E E there is a unique y E F such that x — y E 2ttZ. This is the 
translation congruence property of wavelet sets. Also, for (a.e.) x E E there is a 
unique z E F such that - is an integral power of 2. This is the dilation congruence 
property of wavelet sets. (See section 2.5.6.) 

There is a natural closed-form algorithm for the interpolation unitary Vj^^ 

which maps the wavelet basis for ipE to the wavelet basis for ipp. Indeed, using both 
the translation and dilation congruence properties of {E,F}, one can explicitly 
compute a (unique) measure-preserving transformation a := <t^ mapping R onto 
R which has the property that Vjf ^ is identical with the composition operator 
defined by: 

for all / E L2(R). With this formulation, compositions of the maps a between 
different pairs of wavelet sets are not difficult to compute, and thus products of 
the corresponding interpolation unitaries can be computed in terms of them. 

5.1. The Interpolation Map a 

Let E and F be arbitrary wavelet sets. Let a: E ^ F he the 1-1, onto map imple- 
menting the 27r-translation congruence. Since E and F both generated partitions 
of M\{0} under dilation by powers of 2, we may extend cr to a 1-1 map of R onto 
R by defining (t(0) — 0, and 

a{s) = 2"ct(2~"s) for s E 2'^E, uEl. (29) 
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We adopt the notation a^ for this, and caU it the interpolation map for the ordered 
pair iE,F). 

Lemma 5.1. In the above notation, a^ is a measure-preserving transformation 
from M onto K. 

Proof. Let a := a^. Let fi C M be a measurable set. Let il„ ~ ilr\2"E, n 9 Z, and 
let En = 2^"ri„ C E. Then {iln} is a partition of fi, and we have m{(j{En)) = 
m{En) because the restriction of cr to £^ is measure-preserving. So 



m 



n n 

= ^ 2"m(^(i?„)) = J2 2"^(^") 

n n 

= Y^ m(2"S„) = ^ m(r2„) = m{n). D 

n n 

A function /: M ^ M is called 2-homogeneous if /(2s) = 2/(s) for all s G M. 
Equivalently, / is 2-homogeneous iff /(2"s) = 2"/(s), s e K, n e Z. Such a 
function is completely determined by its values on any subset of M which generates 
a partition of ]R\{0} by 2-dilation. So a^ is the (unique) 2-homogeneous extension 
of the 27r-transition congruence E ^ F. The set of all 2-homogeneous measure- 
preserving transformations of R clearly forms a group under composition. Also, 
the composition of a 2-dilation-periodic function / with a 2-homogencous function 
g is (in either order) 2-dilation periodic. We have f{g{2s)) = f{2g{s)) = f{g{s)) 
and g(/(2s)) = g{f{s)). These facts will be useful. 

5.1.1. An Algorithm For The Interpolation Unitary. Now let 

uE--^u,-, (30) 

where if a is any measure-preserving transformation of M then Ua denotes the 
composition operator defined by Uaf = / o cr~^, / G L^(R). Clearly {o-^)~^ = a§ 
and (C/|)* = t/|. We have U^^^e = V'f since a^{E) = F. That is, 

Ue-$e ^^eoctf ^ ~1k=^e ° erf = —f^Xp = i^F- 
x/ 2tt x/ Ztt 



Proposition 5.2. Let E and F be arbitrary wavelet sets. Then Ue G C2 {D,T). 
Hence T^^Ue^ is the interpolation unitary for the ordered pair {ipE, ^f)- 



Proof. Write a — a^ and Ua = Ue. We have U^i^E = i^F since (j{E) — F. We 
must show 



U^D'^TiPe = D^TU^iPe, nJeZ. 
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We have 

= C/^2-Te-''2-".^^(2-"s) 

This last term is nonzero iff 2^"s G ^, in which case ct^^(2^"s) = cr|;(2^"s) 
= 2^"s + 27rfc for some fc G Z since cr|; is a 27r-translation-congruence on F. It 
foUows that e^*'°" ^^ ''^ = e~*'^ '*. Hence we have 

{U„Wf^i,E){s) = 2-te-''""'"''V^F(2-"s) 

We have shown C/J" € Ct {D,T). Since UEi^E — V-'f, the uniqueness part of Propo- 
sition 3.1 shows that T~^U§T must be the interpolation unitary for {ipE, ipp)- O 

5.2. The Interpolation Unitciry NornicJizes The Commutant 

Proposition 5.3. Let E and F he arbitrary wavelet sets. Then the interpolation 

unitary for the ordered pair (ipE, V'-f) normalizes {D, T}' . 

Proof. By Proposition l5.2l we may work with U^ in the Fourier transform domain. 
By Theorem 6, the generic element of {D, T}' has the form Mh for some 2-dilation- 
periodic hmction h € L°°(R). Write a = aE ^^^ ^cr = Ue- Then 

U-'MhU,=M^oa-^. (31) 

So since the composition of a 2-dilation-periodic function with a 2-homogeneous 
function is 2-dilation-periodic, the proof is complete. D 

5.2.1. C^{D,T) is Nonabelian. It can also be shown ([5, Theorem 5.2 (iii)]) that 
if E, F are wavelet sets with E ^ F then Ue is not contained in the double 
commutant {D^T}". So since U^ and {D,TY are both contained in the local 
commutant oiUjj f at V'Ej this proves that Cg {D,T) is nonabelian. In fact (see 
[5, Proposition 1.8]) this can be used to show that C^,{D, T) is nonabelian for every 
wavelet ip. We suspected this, but we could not prove it until we discovered the 
"right" way of doing the needed computation using s-elementary wavelets. 

The above shows that a pair {E, F) of wavelets sets (or, rather, their cor- 
responding s-elementary wavelets) admits operator-theoretic interpolation if and 
only if GroupjL/^} is contained in the local commutant Cr {D,T), since the re- 
quirement that U^ normalizes {D, T}' is automatically satisfied. It is easy to see 
that this is equivalent to the condition that for each n G Z, a" is a 27r-congruence 
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of E in the sense that (cr"(s) — s)/2tt E Z for all s E E, which in turn implies that 
a'"-{E) is a wavelet set for all n. Here a ~ <t^. This property hold trivially if a is 
involutive (i.e. a^ ~ identity). 

5.2.2. The Coefficient Criterion. In cases where "torsion" is present, so (c^)'^ is 
the identity map for some finite integer fc, the von Neumann algebra generated by 
{D,r}' and U :— U§ has the simple form 

\ Yl ^^h„U'': hn e L°°(R) with K{2s) ^ h„{s), s e M i , 



and so each member of this "interpolated" family of wavelets has the form 

k 



1 '" 

= ^/i„(s)x<T"(ij) (32) 



n=0 



for 2-dilation periodic "coefficient" functions {/i„(s)} which satisfy the necessary 
and sufficient condition that the operator 



I]^4„C/" (33) 



n=0 

is unitary. 



k 



A standard computation shows that the map 9 sending ^g M/^^ [/" to the 
k X k function matrix (/ly) given by 

hij = /la(ij) ° cr"*+^ (34) 

where a{i,j) = {i + I) modulo fc, is a ^-isomorphism. This matricial algebra is 
the cross-product of {D,T}' by the *-automorphism ad{U^) corresponding to 
conjugation with U^ . For instance, if fc = 3 then 6 maps 

Mh,+Mh^Ul:+Mh,{Ul:f 

to 

hi h2 hs ^ 

/13 O (7^^ hi O (T^^ /l2 O (7^^ I . (35) 

^/i2 o (7^"^ /13 o (y^"^ hi o a^"^ 

This shows that '^ijMh^U'^ is a unitary operator iff the scalar matrix (hij){s) 
is unitary for almost all s G M. Unitarity of this matrix-valued function is called 
the Coefficient Criterion in [5], and the functions hi are called the interpolation 
coefficients. This leads to formulas for families of wavelets which are new to wavelet 
theory. 
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5.3. Interpolation Pairs of Wavelet Sets 

For many interesting cases of note, tlic interpolation map a^ will in fact be an 
involution of M (i.e. a o a ~ id, where a :— a^, and where id denotes the identity 
map). So torsion will be present, as in the above section, and it will be present 
in an essentially simple form. The corresponding interpolation unitary will be a 
symmetry in this case (i.e. a selfadjoint unitary operator with square /). 

It is curious to note that verifying a simple operator equation U^ = I directly 
by matricial computation can be extremely difhcult. It is much more computation- 
ally feasible to verify an equation such as this by pointwise (a.e.) verifying explicitly 
the relation aoa = id for the interpolation map. In [5] we gave a number of exam- 
ples of interpolation pairs of wavelet sets. We give below a collection of examples 
that has not been previously published: Every pair sets from the Journc family is 
an interpolation pair. 

5.4. Journe Family Interpolation Pairs 

Consider the parameterized path of generalized Journe wavelet sets given in Sec- 
tion 4.8 Item (iii). We have 

47r' 



Jb 



32tt 

- — ,-47r-4/3)U 



An , 

-n + (3,-—]U 



Air 

y,7r + /3) U 



Att + 4/3, An + 

7 



where the set of parameters /3 ranges — y < /3 < j. 

Proposition 5.4. Every pair {Jp-^, Jfj^) is an interpolation pair. 

Proof. Let /3i, /32 e [-f, f) with /3i < (32. Write a = a'j"^ . We need to show that 

a^{x)=x (*) 

for all a; e R. Since a is 2-homogeneous, it suffices to verify (*) only for x G Jp^. 
For X G J^j n J/32 we have cr{x) — x, hence a'^{x) — x. So we only need to check 
(*) for X G (J/3i\J/32)- We have 

JfhVfh ^ [-I- + Pi, -I- + P2) U [An + APi,An + AP2). 
It is useful to also write 

Jl32\J0i ^ [-47r + 4/3i, -An + 4/32) U [tt + /3i, tt + /32). 

On [—n + /3i, —n + (32) we have cr(x) ~ x + 2n, which lies in [n + f3i,n + /32). 
If we multiply this by 4, we obtain Aa{x) G [An + A(3i,An + 4/32) C J/j^. And on 
[47r + 4/3i,47r + 4/32) we clearly have cr{x) ~ x — Stt, which lies in [— 47r + 4/3i, —An + 
4/32). 

So for X G [— TT + /3i, — TT + (32) we have 

a'^{x) = cr{a{x)) = -a{Aa{x)) = -[4cr(x) - 8n] = (j{x) - 2n = x + 2n - 2n = x. 



Unitary systems and wavelet sets 25 

On [47r + 4/3i,47r+4/32) we have cr(x) = x-Stt, which lies in [-47r+4/3i, -47r + 
4/32). So la{x) e [-7r + /3i,-7r + /32). Hence 



and thus 



<j^x)^4a(^a{x)] =4 



-cr(a:;) + 27r 
4 



= (j{x) + Stt = a; — Stt + Stt = X 



as required. 

We have shown that for all x E Jp-^ we have cr^(x) = x. This proves that 
( J^j , J^jj) is an interpolation pair. 

D 



6. Some Open Problems 

We will discuss four problems on wavelets that we have investigated from an 
operator-theoretic point of view over the past ten years, together with some re- 
lated problems. The set of orthonormal dyadic one-dimensional wavelets is a set of 
vectors in the unit sphere of a Hilbert space H = L^(IR). It is natural to ask what 
are the topological properties of yV(D, T) as a subset of the metric space HI This 
type of question is interesting from a pure mathematical point of view, speaking as 
an operator theorist, and it just may have some practical consequences depending 
on the nature and the degree of depth of solutions. 

Problem 1: [Connectedness] This was the first global problem in wavelets that 
we considered from an operator-theoretic point of view. In [5] we posed a number 
of open problems in the context of the memoir. The first problem we discussed was 
Problem A (in [5]) which conjectured that VV(I?, T) is norm-arcwisc-conncctcd. It 
turned out that this conjecture was also formulated independently by Guido Weiss 
and his group (see [9], [17, [26]) from a harmonic analysis point of view (our point 
of view was purely functional analysis), and this problem (and related problems) 
was the primary stimulation for the creation of the WUTAM CONSORTIUM ~ a 
team of 14 researchers based at Washington University and Texas A&M University. 
(See [26], for the first publication of this group.) This connectedness conjecture 
was answered yes in [26] for the special case of the family of dyadic orthonormal 
MRA wavelets in L^(R), but still remains open for the family of arbitrary dyadic 
orthonormal wavelets in L^(R), as well as for the family of orthonormal wavelets 
for any fixed n and any dilation matrix in R". A natural related problem which also 
remains open, is whether the set of Riesz wavelets is connected. An intermediate 
problem, which is also open, asks whether given two orthonormal wavelets is there 
a continuous path connecting them consisting of Riesz wavelets? (Some evidence 
for a positive answer to this problem is given by Proposition 6.1 below, which 
easily shows that every point on the convex path connecting two wavelets (i.e. 
(1 — t)^ + trj)[s a Riesz wavelet except for perhaps the midpoint corresponding 
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to t = 0.5. Thus this problem has an easy positive solution for many pairs of 
orthonormal wavelets, but no way has been found to get around the midpoint 
obstruction to show that all pairs arc connected, perhaps by some exotic type of 
path.) A subproblem is the same problem but for the set oi frame-wavelets J^{D, T) 
(now widely called framelets). Is the set of all frame- wavelets connected?; or more 
specifically-is the set of all Parseval frame wavelets connected? The reader can 
easily deduce some frame versions of Proposition 6.1 using elementary spectral 
theory of operators which provides some quick-and-easy partial results on paths 
of frames. These are tantalizing, but the main problems still remain open. 

All of these connectivity problems have counterparts for other unitary sys- 
tems. For wavelet systems, they remain open (to our knowledge) for all dimensions 
n and all expansive matrix dilation factors. And for other systems, in particular in 
[14], we showed that for a fixed choice of modulation and translation parameters 
(necessarily, of course, with product < 1) the set of Weyl-Heisenberg (or Gabor) 
frames is connected in this sense, and also it is norm-dense in L^(]R) in the sense 
of Problem 2 below. Although the Weil-Heisenberg (aka Gabor) unitary systems 
are of a simpler operator-theoretic structure than the wavelet systems, and this 
permits the use of some techniques which do not work so well in the wavelet theory, 
even so this perhaps is another reason to think that general connectedness results 
are possible within the wavelet theory. 

Problem 2: [Density] Is the set TZW(D, T) of all Riesz wavelets dense in the 
norm topology in the Hilbert space L^(R)? This was posed as a conjecture by 
Larson in a talk in August 1996 in a NATO conference held in Samos, Greece. 
It was posed in the same spirit as the connectivity problem above, in the sense 
that it asks about the topological nature of TZW{D,T) as a subset of the metric 
space i^(R). Like the connectivity problem it is a global type of problem. A pos- 
itive answer might be useful for applications if it could be given a some type of 
quantitative interpretation. A subproblem of this, which was discussed in several 
subsequent talks, is the same density problem but for the set of frame- wavelets 
T{D^T). (Of course a positive answer for wavelets would imply it for framelets.) 
Like the connectivity problem, this problem makes sense for the family of Riesz 
wavelets for any fixed n and any dilation matrix in R". 

One of the reasons for thinking that this conjecture may be positive is the 
following result, which we think is the most elementary application of operator- 
theoretic interpolation. It is abstracted from Chapter 1 of [5], although the form 
in [5] is a bit different. 

Proposition 6.1. Let hi he a unitary system on a Hilbert space H. If ipi and 'ip2 are 
in W{U), then 

ipi + X'ip2 e nw{u) 

for all complex scalars A with \X\ ^ 1. More generally, ifipi and 'ip2 o-f^ ii^ TZ'W{IA) 
then there are positive constants b > a > such that ipi + A'02 G TZWilA) for all 
A G C with either |A| < a or with |A| > b. 
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Proof. If "01, i/'2 G yV{U), let V be the unique unitary in C^^i^) given by Proposi- 
tion 3.1 such that V1JJ2 ~ ipi- Then 

■V;! +A-02 = (F + A/)-02- 

Since V is unitary, (F + XI) is an invertible element of C^^{U) if |A| 7^ 1, so the 
first conclusion follows from Proposition 3.2. Now assume V'i:V'2 G TZWiU). Let 
A be the unique invertible element of C^^(U) such that Aip2 = "tpi, and write 
■01 + AV'2 = (A + \I)ip2- Since A is bounded and invertible there are 6 > a > 
such that 

a{A) C {z e C : a < |z| < 6} 

where a (A) denotes the spectrum of A , and the same argument applies. 

D 

The above proposition indicates that Riesz wavelets are plentiful. As men- 
tioned above, by writing (1 — t)^p2 + tV'i = ((1 — t)V + *^)V'i) and using the fact 
that the local commutant is a linear space so contains {I — t)V + tl, it follows 
that a convex combination of orthonormal wavelets is a Riesz wavelet except pos- 
sibly for the mid-point corresponding to t ~ 0.5 . So, if ip and 77 are orthonormal 
wavelets, the line in the vector space L^(R) containing the pair tp, rj is in the norm 
closure of the set of Riesz wavelets. Operator theoretic interpolation shows that 
more general linear combinations of finite families of Riesz wavelets are very often 
Riesz wavelets. And if one considers a finite family of wavelet sets C, with union 
S, then the restricted sets of Riesz wavelets, or orthonormal wavelets, or Parse- 
val framelets, or framclets, which are restricted in the sense that they have their 
frequency support contained in S, is always connected if the family of wavelet 
sets is an interpolation family, and these restricted sets of wavelets are very often 
connected even if the interpolation family criterion fails for C. (In fact, it is a 
conjecture that these restricted sets of wavelets are always connected.) Moreover, 
these restricted sets are dense in L'^{S) (considered as a subspace of L^(M)) if the 
family C is an interpolation family, and it is yet another conjecture that they are 
always dense in L^(S). It is also known (see [5] for instance) that the linear span 
W{D,T) is dense in L^(]R). So these facts together suggest that it is probably true 
that the set of Riesz wavelets TZW{D, T) is dense in L^(R). However, the problem 
remains open. Also, as mentioned above under the connectivity problem, one can 
easily deduce some frame versions of Proposition 6.1 using elementary spectral 
theory of operators which provides some quick-and-easy partial results on density 
of Riesz wavelets and framelets. However, the general problem for Riesz wavelets 
remains open. (For the density problem for framelets, we mention that Marcin 
Bownick has recently obtained a significant positive result! It appears that he has 
solved the problem positively for framelets. But apparently it remains open for 
Riesz wavelets.) 

As with the connectedness problem, all of these problems have counterparts 
for other unitary systems. For wavelet systems, they remain open (to our knowl- 
edge) for all dimensions n and all expansive matrix dilation factors (except for the 
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recent interesting franielet density result of Bownick mentioned above). And for 
other systems, as mentioned in the context of Problem 1, in [14] we showed that 
for a fixed choice of modulation and translation parameters with product < 1 the 
set of Weyl-Heisenberg (or Gabor) frames is dense in L^{W) in this sense. This is 
another reason to think that general density results might be possible within the 
wavelet theory. 

Problem 3: [Frequency Support] (See section 4.4.) [Must the support of the 
Fourier transform of a wavelet contain a wavelet set?] This conjecture was posed 
about 10 years ago, by Larson, and the problem still remains open for the case of 
dimension 1 and dilation factor 2. It makes sense for any finite dimension n and 
any matrix dilation, and it apparently remains unsolved in any case. It has been 
studied by several researchers, and Z. Rzeszotnik, in particular, has made some 
progress on the problem. A related problem, also posed by Larson, (see section 
4.4) asks whether a wavelet which has minimal support in the frequency domain 
is necessarily an MSF wavelet. (In other words: Is a minimal support set in the 
frequency domain necessarily a wavelet set). 

Problem 4: [Normalization] (See section 3.3.1.) If {ip,ri} is a pair of dyadic 
orthonormal wavelets, does the interpolation unitary V^ normalize {D,T}'7 As 
mentioned above, the answer is yes if "0 and rj are s-elementary wavelets. This 
problem makes sense for orthonormal wavelets in higher dimensions for matrix 
dilation factors, and for other scalar dilations in one dimension. We know of no 
counterexample for any of these cases. However, the problem might just lie in 
the fact that in most cases, other than for wavelet sets, we have no reasonable 
techniques for doing the computations. This problem was also discussed in context 
in Section 3.3.1, and it could be the most important problem remaining in the 
direction of further development of the unitary system approach to wavelet theory. 
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